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1. Introduction 

The objective of image filtering is to process the 
image so that the result is more suitable then the 
original image for a specific applications. A large 
variety of image processing task can be 
implemented using various filters. A filter that 
attenuates high frequencies while passing low 
frequencies is called low pass filter. Low pass 
filter are usually used for smoothing while high 
pass filter is use for sharpening. 
In literature, Zwirn and Akselrod [1] proposed an 
enhancement method for echocardiography 
images. Echocardiography is such images which 
suffer from noise and difficulty in interpretation. 
Fanga D. and his associated workers [2] proposed 
an image enhancement method for smoothing 
details and preserving edges of an image based on 
nonlinear diffusion equation. Kim B., Kim H. and 
Park D [3] proposed enhancement method for 
fingerprint images. Haddad and Akansu [4] 
employed binomial filters to achieve fast near- 
Gaussian filtering. Clustering filter, an approach 
proposed by Wong [5] which covered the studies 
reveal the use of frequency domain filters. Yusuf, 
Nijad, Sara Tedmory [6] proposed and exploited 
hybrid method for enhancing digital x-ray images. 
In this paper, we performed image sharpening on 
CT and MRI images by Gaussian and Butterworth 
high pass filter and evaluate their performance 
qualitatively (visual inspection) and quantatively 
by calculating mean square error (MSE). MSE 



here is an indication of amount of change in 
original & sharpened images. 

2. Gaussian and Butterworth High 
Pass Filter 

Low pass filters and high pass filters are used to 
show some details in the image while hiding other 
details. Low pass filters blur the image which 
leads to noise reduction while high pass filter 
sharpen some image details, such as edges [6]. In 
the case of Gaussian filtering, the frequency 
coefficients are not cut abruptly, but smoother cut 
of process is used instead. Thus also take 
advantage of the fact that the DET of a Gaussian 
function is a also a Gaussian function. The 
Gaussian high pass filters attenuates frequency 
components that are near to image centre. The 
Gaussian high pass is given by [7] 

Another version of smoothing /sharpening filters is 
Butterworth filter. An advantage with Butterworth 
filter is that we can control the sharpness of the 
filter with the order. A Butterworth high pass filter 
keeps frequencies outside radius Do and discards 
value inside. It has gradual transition from 0 to 1 
to reduce ringing artifacts. A Butterworth filter 
(BHPE) of order n and cut of frequency Do is 
defined as [7] 
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3. Mean Square Error (MSE) 

The Mean Square Error (MSE) is the error metric 
used to compare image quahty. The MSE 
represents the cumulative squared error between 
the compressed and the original image, Mean 
Square Error (MSE) is given by [8, 9] 



MSE = 



I.M,NUi(.'^'n)- 12(771,71) f 



(3) 



M and N are the number of rows and columns in 
the input images. 

4. Methodology 

Image filtering refers to a process that removes the 
noise, improves the digital image for varied 
application. The basic steps in frequency domain 
filtering are shown in Fig. 1. Fourier transform 
will reflect the frequencies of periodic parts of the 
image. By applying the inverse Fourier transform 
the undesired or unwanted frequencies can be 
removed and this is called masking or filtering. A 
filter is a matrix, and components of the filters 
usually vary from 0 to 1. If the component is 1, 
then the frequency is allowed to pass, if the 
component is 0 the frequency is tossed out. 
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Fig.l Basic steps for filtering in frequency domain 

Filtering in frequency domain consist of following 
steps [6] 

1. Pre processing - multiplying the input image 
by (-7 )x+y to centre the transform. 

2. Computing F(u,v), using any of the 
transformation methods. 

3 . Multiplying F(u,v) by 3. filter function H( u, v ). 

4. Computing the inverse transformation of the 
result. 

5. Post processing - multiplying the result by (- 

l)x+y' 




Fig. 2 (a) Original CT image (b) FFT of CT image (c) 
4'^ order BHPF of FFT image (d) GHPF of FFT image 



o 




(a) 




(c) 

Fig. 3 (a) CT image (b) FFT of CT image (c) 4''' order 
BHPF of FFT image (d) GHPF of FFT image 

5. Result and Discussion 

High pass filter give emphasis to the higher 
frequencies in the image. The high pass image is 
then added to the original image so as to obtain a 
sharper image. It may be interesting to experiment 
with width and frequency threshold of the 
Butterworth or the Gaussian high pass filters and 
it will really interesting to compare the sharpening 
which is done in frequency domain and can 
compare it with the sharpening done in spatial 
domain. We will only demonstrate the image 
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sharpening using Gaussian and Butterworth high 
pass filter taking Do=100,n=4(where Do is cutoff 
frequency, n is the order of the Butterworth filter). 
Fig. 2 (a) is the CT image, (b) depicts the FFT of 
the image, and (c) shows the Butterworth high 
pass filter of FFT image (d) shows the Gaussian 
high pass filter of FFT image. Similar examples 
are shown with other CT images as depicted in 
Fig. (3, 4, 5, 6). Equivalent examples can be 
repeated with MRI datasets shown in Fig. 
(7, 8, 9, 10, 11) Now the resultant sharpened 
images of CT and MRI datasets are shown in Fig. 
(12, 13, 14, 15). Now these sharpened images can 
be used in various image processing tasks, like 
edge detection and ridge detection. For 
quantitative analysis, computed the MSE between 
original input images and reconstructed images 
using equation 3. The following calculated values 
for MSE between input images and reconstructed 
images are shown in Table l(a & b) 
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Fig. 4 (a) CT image (b) FFT of CT image (c) 4" order 
BHPF of FFT image (d) GHPF of FFT image 





Fig. 5 (a) CT image (b) FFT of CT image (c) 4" order 
BHPF of FFT image (d) GHPF of FFT image 
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(c) 
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Fig. 6 (a) CT image (b) FFT of CT image (c) 4* order 
BHPF of FFT image (d) GHPF of FFT image 
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(c) (d) (c) (d) 

Fig. 7 (a) Original MRI image (b) FFT of MRI image 

(c) 4th order BHPF of FFT MRI Image (d) GHPF of Fig. 9 (a) Original MRI image (b) FFT of MRI image 
FFT MRI image (c) 4th order BHPF of FFT MRI Image (d) GHPF of 

FFT MRI image 




Fig. 8 (a) Original MRI image (b) FFT of MRI image Fig. 10 (a) Original MRI image (b) FFT of MRI image 
(c) 4th order BHPF of FFT MRI Image (d) GHPF of (c) 4th order BHPF of FFT MRI Image (d) GHPF of 

FFT MRI image FFT MRI image 
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Fig. 14 Sharpened MR images by GHPF 
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Fig. 15 Sharpened MR images by 4th order BHPF 

Table 1 (a) MSB for CT images & (b) MSB for MR 
Images 



MSE Between 
Original CT 
& Sharpened 
CT images 


GHPF 


BHPF 


1 


52.6933 


77.0619 


2 


25.6365 


32.7177 


3 


44.4186 


58.8527 


4 


32.5816 


46.2654 


5 


35.4187 


52.2763 


(a) 


MSE Between 
Original MR 
& Sharpened 
MR images 


GHPF 


BHPF 


1 


87.7435 


117.5782 


2 


29.7703 


46.1320 


3 


27.5597 


45.9106 


4 


18.3468 


29.5037 


5 


60.6688 


88.8303 



(b) 



6. Conclusion 

Gaussian and 4* order Butterworth high pass filter 
with Do=100 is used. High pass filter give 



emphasis on the high frequencies in the image. 
The difference between Butterworth and Gaussian 
filters is that the former provides much sharper 
image than latter. The resultant images by BHPF 
is much sharper than GHPF ,while analysis the 
FFT of CT and MRI images, one sharp spike is 
concentrated in the middle. By applying BHPF & 
GHPF on the images, it is found that the color 
intensities of high pass filter of FFT images 
drastically change. The split edges are elongated 
and wider in BHPF case than GHPF. But in GHPF 
case the split edges are sharper. These elongated 
split edges in BHPF leads to sharper image than 
GHPF. Computing MSE between input images & 
reconstructed images depicts that reconstructed 
images by BHPF possess better image quality than 
GHPF as MSE here indicates change and this 
change is greater in BHPF as compare to GHPF. 
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